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16 Abstract 

The generation of a charge-exchange plasma by a thruster, the transport of this plasma to 
the solar array, and the interaction of the solar array with the plasma after it arrives are all 
described The generation of this plasma can be described accurately from thruster geometry 
and operating conditions. The transport of the charge-exchange plasma was studied experimentally 
with a 15 cm thruster. A model was developed for simple thruster-array configurations. A 
variety of experiments were surveyed for the interaction of the plasma at the solar array. 
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Introduction 


I , 

The use of high-voltage solar arrays can greatly 
reduce or eliminate power processing requirements in space 
electric propulsion systems. The positive high voltage used 
to accelerate beam ions and the low voltage used for the main 
discharge are most promising for direct use of solar-array 
power -- because these two uses represent the largest blocks 
of power in an ion thruster. But both of these uses also 
require substantial areas of solar array to be at high posi- 
tive potential relative to both space and most of the space- 
craft. Such positive potential surfaces, if left exposed, 
could draw excessive electron currents under some conditions. 

The transport of electrons from the ion beam to a 

positive solar array surface was treated by Knauer, et al . 

as an electron space-charge-flow problem.^ Experimental 

measurements were made by Worlock, et al. with an ATS-6 

2 

spacecraft in a large vacuum facility. The ATS-6 spacecraft 
was biased +15 volts relative to the electric thruster neu- 
tralizer, which resulted in substantial electron currents to 
the spacecraft. The current measurements of Worlock, et al. 
were much higher than could be explained by the model of 
Knauer, et al., which indicated the charge-exchange plasma 
surrounding the ion thruster was serving as a conducting 
medium. The ions of this charge- exchange plasma are produced 
by beam ions passing near escaping neutrals, and the transfer 
of electrons from the neutrals to the ions. 
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The transport of electrons from the neutralizer to nearby 
positive surfaces is treated herein as resulting from the in- 
termediate charge-exchange plasma. This electron transport 
can therefore be analyzed in terms of: (1) the production of 

charge-exchange plasma by the thruster, (2) the transport of 
this plasma from the thruster to the solar array, and (3) 
the interaction of the solar array with the plasma after it 
arrives at the array. 
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II. Generation o£ Charge -Exchange Plasma 


The generation of a charge-exchange plasma by an ion 
thruster was described in the preceding annual report of 

3 

this Grant. This process is also summarized herein for 
completeness . 

The overall production rate of charge-exchange ions 

4 

within the ion beam was first calculated by Staggs, et al. 

As will be shown in the next section, electric fields within 
the ion beam are such that almost all charge-exchange ions 
generated downstream of the accelerator grid will escape to 
form the surrounding charge-exchange plasma. 

The following calculation of charge-exchange ion pro- 
duction is similar to that used by Staggs, et al. but is 
presented in a more general form. The neutrals leaving the 
accelerator system are in free molecular flow. The distri- 
bution of neutrals can thus be closely approximated by the 
flow of the same amount of neutral propellant through a sharp 
edged orifice with a diameter equal to beam diameter. Ions i 
pass near the largest number of neutrals by leaving along 
the axis of this orifice. For a conservative assumption, 
then, all ions are assumed to leave on the orifice (beam) 
axis. The integration of neutral density over distance 
along this axis yields 


CO CO 
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where x is the distance downstream of the orifice, r^ is the 

radius of the orifice for beam) , and n is the reservoir 

o , r 

density upstream of the orifice. This density n^ ^ is a 
calculated value that gives the correct loss rate of neutrals. 



N 

o 


n V /4 
o , r o' 


( 2 ) 


where v is the average neutral velocity /8kT /urn . The 
o ° ■' o o 

charge-exchange production rate is thus 


N 


ce 


n r, a N . / 2 
o,r b ce i 


(3) 


with the charge-exchange cross section, 
terms of ion beam current and propellant 
the last equation becomes 


Expressed in 

utilization n , 
u ’ 


N 


ce 


2 

TT r, n q 
b u^ 



(4) 


with q the magnitude of electronic charge. With numerical 
constants substituted, including a typical value of 500°K for 
T^ , the production rate of charge-exchange ions can be written 


N =7.62 
ce 


X 10^^ Jt,^(l-n )a *^r,n 
b ^ 0 “^ ce b u 


(5) 


where beam current J 

. . 2 
^ IS in m 
ce 


b 


is in amperes, charge-exchange cross 
atomic weight A is in atomic mass units. 


section a 



5 


and beam diameter is in m. Typical propellants are mercury, 

cesium, xenon, and argon, for which the charge -exchange cross 

sections at 1000 eV are about 6 x 10 ^ ^ 2 x 10 

4.5 X and 2.5 x 10 ‘ ^^m^ . ^ ^ ® These values change 

- 19 

slowly with energy, with mercury increasing to only 8 x lO 
2 

m at 100 eV. Substituting the charge-exchange cross section 
and atomic weight for mercury gives a production rate of 

N„ = 6-5 . (6) 

Eqs. (5) and (6), then, give the total production rate of 
charge-exchange ions in terms of thruster parameters. As 
these ions escape from the beam volume they carry along 
electrons to form the charge -exchange plasma surrounding the 
ion beam. 
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III. Transport of Charge-Exchange Plasma 

Except for the use of thruster baffles, the transport 
of the charge-exchange plasma from the thruster to the solar 
array was also described in the preceding annual report of 

3 

this Grant. This transport process is also summarized 
herein for completeness. 

The electron density within the ion beam obeys the 
"barometric” equation, 

"e = "e,ref Exp[-qV/kT^] , (7) 

which was introduced by Sellen, et al.^® and verified by 
11 12 

Ogawa, et al. ’ The potential V is defined as zero at the 
reference electron density n Ogawa, et al. found the 

electron temperature within the ion beam is, in electron 
volts, equal to about 0.3 of the injection Cor coupling) 
voltage. Because the Debye length is small compared to beam 
dimensions, the electron density is essentially equal to the 
ion density throughout the beam. Using Eq. C7), equipotential 
contours therefore correspond to contours of equal electron 
density. 

The charge-exchange ions are formed at thermal velocities, 

hence are rapidly accelerated by the potential differences 

in the ion beam. Some detailed trajectories of charge-exchange 

13 

ions have been examined by Komatsu, et al., including the 
effect of initial thermal velocity direction. The approach 
taken herein, though, is to study overall plasma properties 
rather than detailed trajectories. 
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The vacuum facility used was the 1.2-m diameter, 4.6-m 
long chamber at Colorado State University. The thruster was 
a 15-cm SERT-II design, except for the use of dished grids 
with a much higher perveance than the original flat grids. 

All thruster operation was with the screen grid at +1000 V 
and the accelerator grid at -500 V. Conventional Langmuir 
probes were used within the ion beam. For the lower density 
regions outside of the ion beam, 15-cm long probes were used. 
This 15-cm length was divided into three 5-cm lengths, with 
the probe data obtained from the center length. The two end 
lengths were operated at the same potential to serve as 
guards. The thick sheath procedure used to reduce the probe 
data was described by Isaacson. This procedure uses the 
theory and methods of Chen for the accelerating field case. 
This theory is two dimensional, which is why a guarded probe 
was used in low density regions. To minimize cooling of 
electrons by the facility walls, a beam target was operated 
at +60 V relative to the facility during most tests. Some 
tests were also made with a grounded target, but the effect 
of this change on charge-exchange plasma density was small. 
Experimental Data 

The basic problem of interaction between the ion thruster 
and a positive-potential surface is indicated by Fig. 1. The 
simple simulated array was an annular ring surrounding the 
thruster (and in the plane of the accelerator system) with 
an inside diameter of 0.26 m and an outside diamter of 0.92 m. 




Figure 1. Effect of Potential on Current Collected by 
Simulated Solar Array. 
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The electron current collected rises rapidly for small 
positive potentials (relative to the target) . Such a close 
proximity to the thruster is obviously a poor location for 
a solar array, but Fig. 1 clearly shows that excessive cur- 
rents will be collected by positive surfaces near the 
thruster. The problem of ion collection is much less serious 
than that of electron collection, as is also shown by Fig. 1 
(note double current scale) . 

There is no sharp dividing line between the ion beam 
plasma and the surrounding charge-exchange ion plasma. A 
separation between the two regions can be made, though, on 
the basis of relative densities for the two species. Fig. 2 
shows how experimental data are used to find the approximate 
boundary between the two plasma regions. This boundary is 
shown in the surveys of Figs. 3 and 4 by dashed lines. The 
similarity between equipotent ial lines and equidensity con- 
tours is clearly shown within the ion beam region near the 
axis, but less clearly in the surrounding region. The density 
and plasma potential decrease in the radial and downstream 
directions within the ion beam. These variations are typical 
for thruster ion beams. This means that charge-exchange 
ions formed within the ion beam will be directed in the 
radial and downstream directions, and only those that are 
formed close to the negative accelerator grid will be col- 
lected thereon. The assumption was made in the previous 
section that all charge -exchange ions produced in the beam 
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Figure 2. Separation o£ Plasma into Ion-Beam and 
Charge- exchange Regions. 
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escape to the surrounding volume. This assumption is 
supported by the potential distribution within the beam. 

The measured electron temperatures averaged about 5 eV 
in the ion beam and roughly half that value in the surround- 
ing charge-exchange plasma. The electron temperature, 
though, was felt to be the least reliable of the plasma 
properties obtained from probe traces. As mentioned at the 
beginning of this section, the plasma potential is related 
to the plasma density by the barometric equation [Eq. (7)] 
and the electron temperature. The potent ial -density plot 
of Fig. 5 indicates a temperature of about 5 eV within the 
ion beam and about 2.5 eV in the surrounding charge-exchange 
plasma. The lower data in Fig. 5 indicate that the baro- 
metric equation also applies to the charge-exchange plasma, 
but at a lower temperature than for the ion-beam plasma. 

This agreement with the barometric equation should only be 
considered approximate for the charge-exchange plasma, be- 
cause two regions at different electron temperatures cannot 
both satisfy the barometric equation and still have continuous 
densities across an extended common boundary. Because the 
ion-beam plasma is at a higher density and has a higher con- 
ductivity, any discrepancy between the two regions would be 
expected to resolve towards the ion-beam plasma which has 
already been found to agree with the barometric equation to 
a high degree of accuracy. 



Plasma Potential, volts 


Figure 5. Check on Barometric Equation. 
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Transport Model 

A conservative model has been obtained for the trans- 
port of the charge-exchange plasma and compared to surveys 
such as shown in Fig. 3. Conservative means that predicted 
electron densities and currents should err on the high side, 
rather than the low. The model is first developed for an 
isotropic case, then the effect of angular dependence is 
considered . 

Isotropic Model . The total production rate for charge- 
exchange ions (presented earlier) is used with a radial 
variation in ion density. The minimum velocity that could 
be expected (for maximum possible ion density) is the Bohm 
minimum value for a stable sheath, 


V 

ce 


At /m. 
e 1 


( 8 ) 


where k is the Boltzmann constant, T^ is the electron tempera- 
ture (in °K) , and m is the mass of the ion (in kg) , which is 
essentially the mass of a propellant atom m^. The density of 
charge - exchange ions at a radius R (in m) is, therefore. 


n ^ = N ^/4 ttR^v 
ce ce ce 


C9) 


The density of electrons equals that of the charge -exchange 
ions. A positive surface will collect all the electrons that 
arrive at that surface. This electron current density is 

j = n V q/4 , 

■’e ce e ^ ’ 


( 10 ) 
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where is the average electron velocity in the charge- 

exchange plasma, v^SkT ' /ffm . Note that T ' is the electron 

© 0 6 

temperature in the charge-exchange plasma, while is that 

in the ion beam. I£ we use the experimental observation that 

T ' = T /2, then v^ can be expressed as 2/kT /irm . With the 
e e e ^ e e 

numerical constants substituted, the charge-exchange ion den- 
sity and electron current density can be written 


n 


1.49 X 10^^ Ji,^(l-n )a A 
b ^ u'^ ce 


ce 


r, n /T T 
b u o e 


( 11 ) 


2.62 X 10^^ Ju^(l-^ ) 0 A 
b ^ ce 


r, R n vf~ 
b u o 


( 12 ) 


With the substitution of 500°K for T^, 5 eV for (58,000 °K) , 
as well as the atomic weight and charge -exchange cross 
section for mercury, the preceding equations become 


n 


ce 


= 3.3 X 10^2 Jb^^^’^u^/"'b ^^'"u ’ 

j = 0.14 J,^(l-n )/r, n 
-'e b ^ u b u 


(13) 

(14) 


It should be pointed out that lower values of electron 

13 

temperature were obtained by Komatsu, et al., but the value 
found here (5 eV) is consistent with the usual mercury neu- 
tralizer injection voltage of about 20 and the ratio of 0.3 
for electron temperature divided by this voltage, which was 
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11-12 

found by Ogawa, et al. An uncertainty in electron 

temperature, though, will have little or no effect on elec- 
tron currents. Note that in the simple model derived herein, 
the effect of on charge-exchange ion density is cancelled 
by the effect of T on electron velocity, leaving no net 

effect of T on electron current density. 

e ' ' 

Angular Dependence Model . The hemisphere upstream of 
the ion-beam direction is of most interest for spacecraft 
interactions. Charge -exchange ions can leave at 90° to the 
ion-beam direction without interference from the thruster 
body, so that the isotropic model will be used for this 
direction. For the bending of charge-exchange trajectories 
in the upstream direction, we use an experimental observation 
Several surveys of the charge- exchange plasma have shown that 
equipotentials near, and upstream of, the accelerator system 
are approximately normal to the beam direction. The electric 
field in this region is thus nearly antiparallel to the ion- 
beam direction. Using the same value of v'kT /m. as the 
initial velocity of charge -exchange ions in the 90° direction 
this electric field will accelerate ions in the upstream 
direction by an additional potential difference of AV. The 
90° and 180° Cto ion-beam direction) velocity components are 
related by 


2qAV/kT^ = v^g^/v^Q = ctn^e , (15) 


for the 90 to 180° range. The barometric relationship gives 
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\e/\e,90 = Exp[-qAV/kT^'] . (16) 

Again, using half T for T ' , we have 

0 0 

"ce/"ce,90 = Exp [ - 2q6V/kT_, ] . (17) 

Substitution of Eq. (15) into Eq. (17) yields 

"ce''"ce,90 ' , (18) ' 

with 6 restricted to 90 to 180° from the ion-beam direction. 
Inasmuch as electron current density depends linearly on 
electron density, the current density should also be obtained 
by multiplying the value at 90° by the ratio n /n 

0 0 0 0 j ^ U 

Comparison with Experiment . The simple transport model 
described above is compared to experimental results in Figs. 

6 and 7. Inasmuch as the charge-exchange ions of interest 
originate downstream of the accelerator grid, the center of 
the coordinate system for the model should also be shifted 
in this direction. As indicated in the captions of Figs. 

6 and 7, a displacement of one beam radius was used. 

The agreement between experiment and theory in Figs. 6 
and 7 is quite good, considering there are no adjustable 
contants in the theory. The theoretical curves fall above 
the experimental data, which is to be expected from the con- 
servative nature of the assumptions. Similar agreement was 
found for surveys other than the one used for Figs. 6 and 7. 



Electron Density , No./cm' 
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Figure 6. Comparison o£ Theoretical and E:^erimental 
Electron Densities Normal to Ion-beam, Axial Position 
is 7.5 cm Downstream of 15 cm Thruster. 
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Figure 7. Comparison of Theoretical and Experimental 
Electron Densities at Radius of 34 cm from Axial Loca- 
tion 7.5 cm Downstream of 15 cm Thruster. 
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Effects of Thruster Baffles 

If charge-exchange ions could be deflected or collected 
near the thruster, the use of such baffles could be an im- 
portant spacecraft design tool. Tests were conducted using 
as ion beam baffles: (1) a solid cone, (2) a single layer 

screen cone, and (3) a double layer screen cone. The screen 
used had rectangular openings of 0.47 by 0.63 mm and an open 
area fraction of 0.59. The intent was to collect charge- 
exchange ions that were formed well downstream of the 
thruster. The use of screens, of course, was to facilitate 
the escape of neutrals, hence reduce the charge-exchange 
production relative to the solid cone. 

The method of selection for the baffle location is 
indicated in Fig. 8. The baffle was placed just outside the 
region dominated by beam ions. The length Z was made 20 cm 
to assure that most of the charge-exchange process took 
place within the cone. For the double-layer screen cone, the 
inner screen was located as shown in Fig. 8. The actual 
baffle construction is indicated in Fig. 9. The double- 
layer screen cone was identical in construction to the single- 
layer cone, except that a second layer of screen surrounded the 
outside of the support members. The gap between the two screen 
layers was 2.5 cm. 

The plasma densities measured for the three different 
baffles are shown in Figs. 10 - 12. As mentioned earlier, 
the target was maintained at +60V relative to the vacuum 
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Figure 10. Map o£ Electron Density for Solid Cone. 



(b) Baffle at +30 V relative to facility 













Figure 11. Map o£ Electron Density for Single-Layer Screen Cone. 
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(b) Baffle at +30 V relative to facility 






Figure 11. Concluded. 
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Figure 12. Map of Electron Density for Double-Layer Screen Cone. 
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(b) Baffle at +30 V relative to facility 
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(c) Baffle at +60 V relative to facility 
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facility to avoid ion reflection at the facility walls. The 
baffles were biased at 0, +30, and +60V relative to the 
facility. It was not practical to operate the baffles at 
potentials more positive than the target due to the large 
electron currents that would be collected. 

The performance of the baffles operated at the dif- 
ferent bias voltages is indicated in Fig. 13 for a single 
location relative to the thruster. It had been hoped that 
baffle potentials negative, relative to the target, would 
tend to collect the charge -exchange ions. As shown in Fig. 
13, though, there was no clear beneficial effect of such a 
bias. Performance was also evaluated at other locations 
upstream of the thruster, but the results were qualitatively 
similar to those shown in Fig. 13. 

Inasmuch as bias power supply would be an added com- 
plication, the bias supply cannot be recommended in the 
absence of a clear advantage for its use. The case of 
baffles at spacecraft potential was approximated at the 
+60V bias (baffles at target potential) . This condition 
was used for comparison of Fig. 14, which was made at a 
radial distance of 30 cm. The use of baffles significantly 
decreased the plasma density upstream of the thruster, 
where most, or all, of the spacecraft would be located. 

The largest decrease was found with the solid cone, while 
the smallest decrease was with the single layer screen cone. 
The poorer performance of the screen cones indicates that 
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Figure 13. Effects of Cone Potential and Construction 
at a Radius of 30 cm and an Axial Location 17.5 cm Up- 
stream of Accelerator Grid. 



Electron Density , No. /cm 



Figure 14. Variation of Charge- exchange Plasma Density at a Radial Distance 
of 30 cm. Axial Distance Increases Downstream and is Measured from Accelerator 
Grid. 
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the escape of neutrals through the screen was more than 
offset by the simultaneous escape of charge-exchange ions. 

The screen holes were smaller than the local Debye length, 
so that a finer mesh size would not be expected to change 
the results. 

At least some of the improvement due to the solid cone 
in Fig. 14 resulted from simply moving the source of the 
charge-exchange ions further downstream. The effect of this 
source displacement is removed in Fig. 15. Zero axial 
distance for the no baffle configuration was taken as one 
beam radius downstream of the accelerator system, as men- 
tioned in connection with Figs. 6 and 7. No rigorous selec- 
tion of a single distance can be made, because the effective 
source location depends on the specific potential distribu- 
tion within the ion beam. One beam radius, however, seemed 
a reasonable compromise for the no baffle conditions in- 
vestigated. In the derivation of the simple model presented 
earlier, the effect of the ion-beam radius is through the 
neutral distribution. With the solid cone, then, the ef- 
fective beam radius should be the exit radius of the cone. 

The zero axial distance for the solid cone was therefore 
taken as one cone exit radius downstream of the cone. 

The near superposition of peaks for the two sets of data 
indicates reasonable assumptions were made for effective 
source locations. 

The solid cone data of Fig. 15 are everywhere above 
the no baffle data. This result indicates that any collection 



Axial Distance , cm 


Figure 15. Variation of Charge- exchange Plasma Density at a Radial Distance 
of 30 cm. Axial Distance Increases Downstream and is Measured from One Beam 
Radius Downstream of the Accelerator Grid for no Cone, and from One Cone Radius 
Downstream of Cone Exit for Solid Cone. 
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of charge-exchange ions is more than offset by the higher 
production caused by higher neutral densities within the cone 
The major effect of a baffle thus appears to be the movement 
of the charge-exchange ion source in the downstream direction 
The improvement due to this effect is partially offset by 
an increase in total charge-exchange ion production. 
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IV. Interaction of Solar Array with Plasma 

A specific spacecraft configuration is useful in evalua- 
ting the interaction of thruster and solar array. The 
configuration assumed herein is shown in Fig. 16. The 
thrusters are assumed to be 30 cm models with a maximum 
beam current of 2 amperes and a propellant utilization of 
0.9. The total array power is 25 kw at one astronomical 
unit from the sun, and a maximum of eight thrusters is 
assumed to operate at any one time. The distance from the 
solar array to the thruster exhaust plane is varied from 1.2 
to 2.4 meters. 

The charge-exchange plasma properties near the solar 
array were calculated using the model presented earlier. 

The point of origin for this plasma was assumed to be one 
beam radius (15 cm) downstream of the center thruster. The 
case of maximum plasma density was obtained with the effect 
of one operating thruster multiplied by eight for the total 
thruster array. This approach has the implicit assumption 
that the thrusters are far apart. For a closer spacing, 
there would be an additional contribution due to thruster 
ions from one operating thruster passing through the neutral 
effluxes of the other operating thrusters. 

Unprotected Solar Array 

The electron density, saturation electron current density, 
and Debye length are shown in Figs. 17-19 as a function of 
location on the solar array and the distance of the array 




Electron Density , No./cm 
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Figure 17. Distribution of Charge -exchange Plasma Density at Solar-array Panel. 
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Figure 18. Distribution of Current Density at Solar-array Panel. 
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Figure 19. Distribution of Debye Length at Solar-array Panel. 
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upstream of the thruster exhaust plane. The current density 

was numerically integrated over the array area and shown 

by the solid lines in Fig. 20 as a fraction of solar array 

current (25 kw divided by array voltage) . This integration 

assumed all of the solar array was positive relative to the 

thruster neutralizer, so that all the electrons arriving at 

the array would be collected. There is also an effect of 

plasma sheath thickness. As the array becomes more positive, 

the sheath thickness increases. Comparison of Child's law 

and Debye length equations will show that the ratio of Child's 

3 / 

law distance to the Debye length is 1.26 x (V /E ) ^ where 

3 . 0 

V is the array voltage relative to the plasma and E is the 
electron temperature in eV. This relationship assumes a 
plane parallel geometry of sheath and array, but should 
approximate the effective increase in solar array width due 
to sheath effects. With the area increase due to sheath 
effects, the integrated electron current is shown as a 
fraction of array current by the dashed lines in Fig. 20. 
Inasmuch as the entire array was assumed to be at maximum 
potential, this sheath thickness correction should be a 
worst case value. 

The ratio of collected electron current to array cur- 
rent does not translate directly into power loss. Because 
the return circuit is through the neutralizer, different 
cells are loaded with different values of parasitic current. 
Although the current- voltage load curve is usually the same 


Collected Electron Current os Fraction 
of Solor Array Current 
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Figure 20. Electron Collection Fraction for Solar Array. 
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for all cells, the varying parasitic currents result in a 
range of operating points on this load curve. The overall 
curve therefore does not shift to lower current values 
but is distorted by the mismatch into a new curve shape. 
Depending on the allowable voltage change for the array, the 
power loss fraction can be substantially less than simply the 
fraction of array current that is collected from the surround- 
ing plasma. 

A number of simplifying assumptions were made for Fig. 

20, but some quantitative conclusions can still be drawn. 

It should be kept in mind, though, that these conclusions 
depend on the validity of the transport model, which has been 
verified over only a small range of radius ratio. From Fig. 

20 it appears that an unprotected array operating at +200 
volts should experience a small ('^'lO percent) loss due to 
interaction with a thruster, or thrusters. With a center tap 
grounded to the spacecraft (neutralizer potential), a 400 volt 
array should have a minimal C'^^5 percent) loss. The interaction 
losses can be further reduced by placing the highest potential 
cells at the tips of the arrays, where the parasitic current 
density is smallest. The penalty due to parasitic electron 
currents increases rapidly with array potential. At 1000 
volts, which is about the value required for acceleration of 
beam ions, the parasitic current could be 20 to 50 percent of 
the array current. Parasitic currents of these magnitudes 
can result in substantial power losses. For more precise 
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evaluations of unprotected arrays, though, detailed cal- 
culations must be made using the actual distributions of 
potential on the arrays. 

More precise evaluations should also include a more 
accurate approach to sheath thickness. In the calculation 
used herein, the one dimensional sheath thickness (Child's 
law distance) was added around the edge of the projected 
solar-array area. Because the added area is found at the 
edge of the solar array, a radial - inf low solution should be 
used in place of the planar solution (Child's law) used 
herein. A radial - inf low solution, however, will depend on the 
specific geometry of the array edge and the surrounding 
plasma sheath. The use of segments of radially symmetric 
solutions should give intermediate accuracy, and hence 
indicate the error associated with the use of a planar 
solution length. 

Although the effects shown in Fig. 20 were calculated 
for a specific spacecraft, they are approximately valid 
for other spacecraft sizes. If all the spacecraft dimensions 
were doubled, for example, the solar power would be in- 
creased by a factor of four. This increase would permit 
four times as many 30-cm thrusters to be operated at one 
time, resulting in four times the maximum generation rate 
for charge-exchange ions. The inverse-square variation 
in charge -exchange plasma density would thus result in the 
same mean values of electron density and saturation electron 
current at the solar array. Assuming the same array voltage. 
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the electron currents collected on the original and scaled- 

up array would be the same fraction of the array current. 

The plasma sheath correction would, of course, fail to scale 

with a size change. (It would give a small performance 

increase for the large array.) 

Protection of Solar Arrays 

An obvious possible solution for the interaction at 

high voltages is to cover the array with a layer of insulation, 

thereby isolating the array from the surrounding charge- 

exchange plasma. A variety of materials -- polyimide (Kapton) , 

flourinated polyethylene (Teflon), and glass -- will withstand 

high voltages in thin layers. A large area of insulation 

should be expected to have a few small holes, though, either 

from defects in manufacturing or collisions with micro- 

meterorites. The effectiveness of insulation for a solar 

array thus reduces to the effects of these small holes. 

From Fig. 17, the range of interest for density is from 

about 10^ to 10^ electrons/cm^ . A number of studies have 

been made of holes in insulators in a plasma environment, 

some of which are in this range of density. The current 

collected through an 0.38 cm hole in an insulator sheet was 

1 7 

evaluated by Kennerud and is shown in Fig. 21. The satura- 

- 7 2 

tion electron current density is only about 0.1 x 10 A/cm 
for the conditions studied by Kennerud, so that the large 
observed currents must come from surface leakage and/or the 
focusing effect of the adjacent plasma sheath. Similar results 
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Figure 21. Leakage Current Collected by an 0.38 mm hole in 0.13 mm Thick 
Polyimide (Kapton). From Ref. 16. 
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4 2 18 

were also shown by Grier and Domitz at 10 electrons/cm . 

At higher densities 10° electrons/cm ) , investigations by 
19 2 0 

Cole, et al. and Grier and McKinzie showed high currents 
accompanied by hole damage and enlargement. 

Small holes are thus found to be very effective in 
collecting electrons from the charge-exchange plasma. The 
effectiveness is sufficient that the protection by a nearly 
continuous layer of insulation is questionable. In fact, 
the concentration of electron currents at just a few holes 
may be far more damaging (as far as electron currents are 
concerned) than the more uniform current distribution of 
an unprotected array. Solar-cell cover glasses, with ex- 
posed connections between cells, are probably equivalent 
(for electron collection) to unprotected arrays. 
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V. Concluding Remarks 

The environmental interaction between electric thrusters 
and high-voltage solar arrays is due to the charge-exchange 
plasma generated by the thrusters, with the most serious 
effect being the conduction of electrons to positive sur- 
faces. The charge-exchange ion production can be readily 
calculated from thruster dimensions and operating parameters. 
A simple model is given for the transport of these ions from 
the thruster, or thrusters, to the solar array. Attempts 
to deflect or collect these charge -exchange ions at the 
thruster were largely ineffective in reducing the plasma 
density upstream of the thruster. The only significant 
effects of these attempts was the displacement of the 
effective source of charge -exchange plasma in the downstream 
direction. Attempts to protect a solar array with a layer 
of insulation also appears ineffective, with large currents 
observed through small holes. 

The investigations to date have been preliminary in 
the sense that parameters have not been varied in a thorough 
and systematic manner. Further study would therefore be 
expected to improve both the accuracy and completeness of 
thruster/array interactions. For example, the transport 
model presented was verified over only a limited radius 
ratio. Investigation over a wider range of this parameter 
may result in a new transport model. New concepts may also 
significantly alter the interaction problem. As an example 
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of a possible new concept, a conducting grid, printed on 
the plasma side of an insulating layer, should break up 
either surface conduction or plasma focusing to small holes. 
The current to a hole should therefore correspond to 
approximately the mesh size of the conducting grid. The 
sizing of the grid would, of course, depend on both the 
performance loss due to the presence of the grid and the 
probability of holes - due to either manufacturing defects 
and micrometeorites. Because of the parametric and 
concept limitations, the study presented herein should be 
considered only a preliminary description of the thruster- 
array interaction. 
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